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ABSTRACT 

In the last few years, over 43 millisecond radio pulsars have been discovered by targeted searches of uniden- 
tified 7-ray sources found by the Fermi Gamma-Ray Space Telescope. A large fraction of these millisecond 
pulsars are in compact binaries with low-mass companions. These systems often show eclipses of the pulsar 
signal and are commonly known as black widows and redbacks because the pulsar is gradually destroying its 
companion. In this paper, we report on the optical discovery of four strongly irradiated millisecond pulsar com- 
panions. All four sources show modulations of their color and luminosity at the known orbital periods from 
radio timing. Light curve modelling of our exploratory data shows that the equilibrium temperature reached 
on the companion's dayside with respect to their nightside is consistent with about 10-30% of the available 
spin-down energy from the pulsar being reprocessed to increase the companion's dayside temperature. This 
value compares well with the range observed in other irradiated pulsar binaries and offers insights about the 
energetics of the pulsar wind and the production of 7-ray emission. In addition, this provides a simple way of 
estimating the brightness of irradiated pulsar companions given the pulsar spin-down luminosity. Our analy- 
sis also suggests that two of the four new irradiated pulsar companions are only partially filling their Roche 
lobe. Some of these sources are relatively bright and represent good targets for spectroscopic follow-up. These 
measurements could enable, among other things, mass determination of the neutron stars in these systems. 
Subject headings: binaries: general — pulsars: general — pulsars: individual (PSRs J1810H-1744, J0023H-0923, 
J2215H-5135, J2256-1024, B1920H-57, J1023H-0038) 



1. INTRODUCTION 

The Fermi Gamma-Ray Space Telescope (hereafter simply 
Fermi) has discovered hundreds of unidentified point sources 
whose positional uncertainties are small enough to enable 
deep targeted pulsar searches with the world's largest radio 
telescopes. These searches have resulted in the discovery 
of 43 energetic millisecond pulsars (MSPs) to date (|Cognard 



et al. 201 1 ; Hessels e t al.|201 1 [ 'Kerr et a l.|2012t [Ransom et aL 
201 1 , Ray et al. 2012!). Intriguingly, many of these are "black 
widow" and "redback" systems, of which previously only a 
handful were known - Fermi has so far increased their num- 



ber to over twenty ( Roberts 
Black widows are energetic 



SPs (spin-down luminos- 



ity E ^ 10 erg s ) with very low mass (few 0.01 M©) 
companions in compact, few-hour orbits. The original black 
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widow, PSR B1957H-20 ( jFruchter et al.[[1988b| ), and other 
members of the class derive their name from the fact that 
their strong relativistic wind ablates the surface of the com- 
panion star, whi ch might have bee n significantly more mas- 
sive in the past ( Fruchter et al.|1 98^b ). Redbacks, on the other 
hand, have similar pulsars but more massive companions, of 
^O.2M0. As yet, is it not clear whether redbacks are the 
close progeny of accreting X-ray systems, nor whether they 
can evolve into black widow systems. The evidence that the 
prototype redback PSR J1023-I-0038 accreted matter no t long 
ago is suggestive of the former at least (Ar chibald et al.|2 009 ). 

Because of their "cannibalistic" behavior, it has been pro- 
posed that some blac k widows eventually completely ablate 
their compa nion ( Klu zniak et al.|198"8 van den Heuvel & van 
Paradijs'1988 ). If this is the case, it could at least partially re 



solve the conundrum of the existence of isolated MSPs, since 
these sources clearly appear to have evolved in binary sys- 
tems despite the fact that they are now isolated. It appears, 
however, that at least for the original blac k widow pulsar, 
PSR B 1957-1-20, the e vaporation is too slow ( jEichler & Levin- 
[son[[1988[|Eichler & Gedalin 1995). Hence one has to won- 
der whether the known black widows are somewhat special 
(and less extreme than the isolated millisecond progenitors) 
or whether the scenario is simply wrong. Recent work sug- 
gested that evolution in a triple system might provide a vi- 
able formation channel to explain a sma ll fraction of the iso- 
lated MSPs and peculiar binary pulsars (Portegies Zwart et al. 
[201 1[ [Freire et al.[[201 1[ ). Nonetheless, the ablation mecha- 
nism still appears as the most plausible general scenario for 
creating isolated MSPs. 

Black widows and redbacks usually display extended 
eclipses at radio wavelengths, which are accompani ed by 
rapid variations of the dispersion measure (Stappers et al. 
[1996| l. At X-ray energies, persistent emission is often vis- 
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ible in addition to pulsations, which likely results from the 



shocked relativistic wind collid ing with the companion ( |Stap 
pers et aL]|2003 Tavani|[r993 1. Black widows and redoacks 



were thought to be a relatively small fraction of the total MSP 
population, but it has become clear that previous wide-field 
pulsar surveys have missed many sources (e.g. because of 
eclipses). In contrast, since energetic MSPs are efficient 7- 
ray emitters, targeted radio searches of 7-ray sources are bi- 
ased towards finding them, especially since repeated searches 
of the same source have a better chance of catching it out of 
eclipse. About one third of the pulsars discovered by target- 
ing 7-ray so urces have been black widows or redbacks (Ray 
et al.||2012 i. In the light of these discoveries, black widows 
might after all offer a viable channel for the formation of at 
least some of the isolated millisecond pulsars. 

Optical and near-infrared observations of black widows and 
redbacks are an important probe of the state of the com- 
panion and the energetics of the system (see, e.g. 



Fruchter 



eraL]|1988a[|van Paradijs"eral . 1988; Stapp ers et al.[|2001| . 
ITie optical light is dominated by the companion - the pul- 



optical light is dominated by the companion - tiie pul- 
sar contributes a negligible amount - and shows significant 
flux and color variations. Mostly, these reflect strong ir- 
radiation of the hemisphere facing the pulsar, but super- 
posed on this are ellipsoidal variations due to the tidal dis- 
tortion. Combined, these allow one to constrain the incli- 



nation and other phy sical parameters ( Stappers et al.| 2001 
[Reynolds et al.l| 2()07| l. Combined with phase-resolved spec 
troscopy, this ca n be used to determine the component masses 
(e.g. Iv an Kerk wijk et al.pOl 1| [Romani & S haw 201 1 ; Ro- 



mani et al. 20lli. Using this technique, van Kerkwijk et al. 



( |2011[) infeEedfliat PSR B1 957-H20 is likely very massive, 
Jm±0.l2MQ. Similai-ly, Romani et alj ( [2012) 1 find evi- 



dence that PSR J131 1-3430 is also heavyweight - perhaps as 
much as ^2.1 Mq. These large masses suggest mass trans- 
fer was relatively effective, in contrast to what is inferred for 
other pulsar binaries (e.g., Lin et al.]|201 1[ Antoniadis et al. 



2012 1, posing both an interestmg quandary as to why this 



might be the case, and an opportunity to probe the upper mass 
limit of neutron stars. 

While this kind of light curve modelling is also po ssible in 
other types o f binary neu tron star sy stems (see, e.g., Munoz-| 



Darias et al.|2009, Wang et al. 201 1[ for low-mass X-ray bina- 
ries), black widows and redbacks are much cleaner systems. 
Indeed, the only source of optical light is the companion since 
there is no accretion disk or je|^ Also, the irradiation is due 
to relativistic photons/particles which penetrate several opti- 
cal depths inside the companion's photosphere. As a result, 
the atmosphere remains in quasi equilibrium and this avoids 
the formation of prominent emission line features. Finally, 
because the neutron star is a pulsar, the radio timing provides 
a set of accurate orbital parameters. 

We have searched for the optical counterpart to four of the 
new black widow/redback systems, using Gemini data com- 
plemented by New Technology Telescope (NTT) and archival 
Swift data. In this paper, we present the optical discovery 
of these companions (Section |2]l and use their light curves to 
constrain the system parameters by modelling them using the 
synthesis code Icarus (Section|3]l. We find that the temper- 
ature increase on the dayside of the irradiated pulsar compan- 
ions typically corresponds to a conversion efficiency ^ 15% of 
the incident energy available from the host pulsars' rotational 

" PSR Jl()23+0038 was observed to have a disk before but it has since 
disappeared (see, e.g., Archibald et al. 2009 1. 



spin-down (Section|4]i. 

2. OBSERVATIONS AND ANALYSIS 

We were awarded 7.67 hr on Gemini North (program GN- 
2010B-Q-77) to search for the counterparts of four ener- 
getic binary MSPs with the GMOS-N instrument ( Hook et al.| 
|2004l i: PSRs J1810H-1744, J0023H-0923, J2215+5135 and 
J2256-1024. The first three pulsars were found in a targeted 
GBT search of unidentified Fermi point sources (Ban gale et 
al. in prep; Hessels et al. in prep; 'Hessels et al.'20ir) while 
the last one was discovered during the 350 MHz GBT pulsar 
drift scan survey (Boyles et al.|20l2| [Lynch et al.|2012( Stairs 
et al. in prep). Table [T] presents the main properties of the 
targets, and a detailed discussion about each source follows in 
Section]?] Our observing program had two main goals: detect 
the pulsar companions and, if they are detected, identify vari- 
ability at the orbital period. Our observing strategy consisted 
of collecting data at four different epochs, with each observ- 
ing session consisting of a 320-s /-band, 620-s g-band and 
320-s /-band exposure sequence. Since we observed under 
non-photometric conditions (cloud cover: 90%, image qual- 
ity: 85%, sky background: 80%, water vapor: any), resolution 
was not an issue and we binned the EEV CCD detector by a 
factor 2 X 2 in order to reduce the readout time to 35 s per 
exposure (fast mode) while still properly sampling the point 
spread function. Using this strategy, we were expecting to 
find inter-epoch variability and, given the short orbital peri- 
ods of these binaries, maybe even intra-epoch variations (for 
the /-band). 

We reduced the data following standard procedures, im- 
plemented using custom Python scripts. We used standard 
Gemini nightly calibration data to remove the bias and flat- 
field our science frames. We registered our frames astromet- 
rically relative to the UCAC3 catalog for PSRs J1810H-1744, 
J0023H-0923 and J2215H-5135, and relative to SDSS for 
PSR J2256-1024 (the only source that faUs within the cov- 
erage of SDSS Data Release 7). For all but PSR J0023H-0923, 
a few tens of reference stars fall within the central frame 
of the CCD detector and hence our calibration yielded posi- 
tional uncertainties dominated by systematic uncertainties, of 
0."07 and 0"1 in UCAC3 and SDSS DR7, respectively. For 
PSR J0023H-0923, only three UCAC3 stai's fall on the central 
CCD chip. In this case, we calibrated the image that looked 
the cleanest using the UCAC3 catalogue, and then calibrated 
the other images in the same band relative to the reference one 
using a list of bright stars found in all images. The measured 
optical positions are reported in Table ]T] and agree with the 
radio position derived from the timing. 

We performed aperture photometry using an extraction ra- 
dius of 5 pixels (i.e. 0.73 arcsec at the plate scale of 0.146 
arcsec per binned pixel), and sky inner and outer annuli of 
10 and 15 pixels, respectively. We calibrated our photometry 
against bright, non-saturated stars appearing in the SDSS DR7 
catalogue for PSR J2256-1024, and USNO-Bl in the case of 
the three other targets. For USNO-Bl, we converted the cat- 
alogue magnitudes from the photographic B, R, and I mag- 
nitudes to th e ugriz system using the transformation of Jordi] 
]eraL] ( ]2006l l, / - / = (0.24 7 ±0.003)(j;- /) + (0.329 ± 0.002), 
for the /-band, and from Lupton] (J2005]), B-g = 0.3130(g- 
r) + 0.227 1 and - r = -0. 1 837Ti^ 7Pa097 1 , for the ^-band. 
In the case of PSR J18 10-1-1744, there is a neighbouring star 
that might contaminate our aperture photometry. For this rea- 
son, we performed point-spread function photometry using a 
Moffat profile f(r) cx (1 + r^/tr^)"'^, with /? = 3 and a opti- 
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Table 1 

Measured and Inferred Source Parameters 



Quantity 


10023+0923'' 


J2256-1024'' 


11810+1744" 


12215+5135" 


B1957+20''' 


11023+0038= 


Measured 


Right Ascension^' 


uu Zj io .oy(,zj 


22'"56'"56^39(1) 


18''10"'37'.28(1) 


ZZ 1 J jZ .Oo(^i ) 


iV jV jO . 1 1 


lU Ld 4/ .Kyi 


Declination'' 


09° 23'24". 18(20) 


-10°24'34".37(12) 17°44'37".38(7) 


51°35'36" .45(10) 


20°48'15".12 


00°38'41".15 


Gal. Longitude (degree) 


111.38 


59.23 


44.64 


99.87 


59.20 


243.49 


Gal. Latitude (degree) . 


-52.85 


-58.29 


16.81 


-4.16 


-4.70 


45.78 


/7,0.1-looGeV 


10.7± 1.5 


10.1± 1.4 


25.5 ±2.1 


10.9± 1.6 


16.7±1.9 


5.4± 1.0 


(IQ-'^ergs"' cm-2) ... 














Pspin (ms) 


3.1 


2.3 


1.7 


2.6 


1.61 


1.69 


Porb (h) 


3.3312 


5.1092 


3.5561 


4.1401 


9.1672 


4.7543 


Tlscnode (MJD (TDB)) . 


55186.11343 


54853.22391 


55130.04813 


55186.16449 


48196.06352 


54801.97065 


;i:(lt-s) 


0.035 


0.083 


0.095 


0.47 


0.089 


0.343 


DM(pccm-3) 


14 


14 


40 


69 


29.12 


14.33 


Av^ 


0.37 


0.14 


0.43 


1.15 






Inferred 


dDM (kpc)f 


0.69 


0.65 


2.00 


3.01 


2.49 


0.62 


M™"(M0)E 


0.017 


0.030 


0.045 


0.213 


0.022 


0.138 


L,j,P (103^ ergs-')'' .... 


1.51 


3.95 


3.97 


5.29 


16.0 


9.82 


Lj (10^3 ergs-')' 


1.28 


1.21 


3.05 


1.30 


2.00 


0.65 


aiReY 


1.27 


1.69 


1.33 


1.53 


2.49 


1.65 


Rl (R©)'' 


0.13 


0.22 


0.19 


0.36 


0.29 


0.34 


Light Curves (minimum / maximum / quadrature)' 


/-band 


24.3/21.7/22.7 


24.3/20.8/22.1 


22.3/ 19.5/20.3 


19.5/18.6/18.9 






g-band 


28.0/23.4/25.0 


28.0/22.2/26.8 


23.2/ 19.2/20.2 


19.9/18.1 / 18.7 






Light Curve Fitting 


inclination (degree) .... 


58 ±14 


68±11 


48±7 


66 ± 16 






filling factor 


0.30 ±0.30"' 


0.40 ±0.20 


0.80 ±0.30 


0.99 ±0.03 






Tfiight (K) 


2900 ±700 


2450 ±350 


~ 4600" 


4800 ±450 


~ 2500 


~5600 


7day(K) 


4800 ± 2000 


4200 ±700 


> 8000" 


6200 ±500 


~ 5800 


~ 6650 


r,„(K) 


4600 


4100 


7800" 


5550 


5750 


5580 


^irr 


0.17 


0.07 


0.60" 


0.15 


0.15 


0.09 



"Timing data (Pspin. ^orbi and L^^ p from Hessels et al. (in prep.). 

''Timing data (Pspin. Rmb. ^. and L^^ p from Stairs et al. (in prep.). 

■^PSRs B 1957+20 and J1023+0038 are shown for comparison. Data from the ATNF Pulsar Catalogue ^Manchester et al.|2005||http : / /www . atnf . csiro . | 
'^u/research/pulsar/psrcat I. 

"Positions derived from optical astrometry. Uncertainties are dominated by the catalog systematics of 0"07 in UCAC3 (J0023+0923, J18 10+1744 and 
J2215+5135) andOf'l in SDSS DR7 (J2256-1024). 
■^Total reddening along the line of sight (see Section|3|. 
Distance based on the dispersion measure NE2001 model JCordes & Lazio : 

eAT- ™ : „^„..™; ^ a \/i i„ n no 1 \ 



2002 



^Minimum companion mass assuming a 1.4Mq pulsar and a yu^ orbital mclinatibn. 
''Spin-down luminosity inferred from Pspj,, and its derivative, and a moment of inertia of 10**' gcm^. 
'7-ray luminosity estimated from the 7-ray flux and dispersion-measure distance. 
J Orbital separation a = a-(1 + 1.4/M™"). 
'^Companion roche radius Rt = 0.46a[M™'"/(M™" + 1 .4)]'/3. 

'The reported magnitudes are inferred from the light curve modelling presented in the text. 

"The probability distribution for this parameter is highly non-Gaussian, hence the value must be taken with caution. The median value is 0.40, the mode at 
~ 0.15 and the distributio n ext ends with a heavy tail all the way to unity. 

"As explained in Section[33] we believe that our analytic estimate is more robust than the numerical values for the case of PSR J1810+17. Hence we present 
the analytic results in this table. 



mized using a set of bright, non-saturated stars for each frame. 
Reference stars were fitted individually while the immediate 
vicinity of PSR J18 10-1-1744 was simultaneously fitted for the 
source and the two other stars located East and North-East 
from it. 

Our complete photometric results are available online in Ta- 
bles l2j [3] |4] and p] Our photometric errors were calculated by 
adding m quadrature the sky background, the photon counting 
noise and the intra-band relative zero-point. The zero-point 
calibration errors correspond to the standard deviation of the 
mean of the zero-point calibration for each band to the refer- 
ence catalog stars, which were added in quadrature to the cat- 
alog systematic calibration. In the case of PSR J0023-I-0923, 
very few catalog stars overlap with our field and hence the 



band calibration is poorer than for the other systems analyzed 
here. The SDSS systematic calibration error is 0.02 mag (for 
PSR J2256-1024), while that of USNO-Bl.O is 0.3 mag (for 
PSRs J1810H-1744, J0023H-0923 and J2215H-5135). 

The field of PSR J22 15-1-5 135 was also serendipitously ob- 
served by the UVOT instrument on-board Swift while it was 
monitoring a nearby gamma-ray burst. We performed photo- 
metric reduction of the publicly available datcpl which were 
all obtained in the uvwl-band (Poole et al. 20(MpAc = 260 nm, 
AA = 69.3 nm). For the UVOT data, we used a 5 pixel 
aperture (i.e. 5arcsec at 1.004 arcsecpixel"^) and sky in- 
ner and outer annuli of 15 and 30 pixels, respectively. We 



''^ http : //heasarc . gsf c . nasa . qo\ 
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took the photometric zero-point from Breeveld et al. ( 201 l| l, 
zPuvwi = 18.95 ±0.03 (AB system), to convert our count rate 
to magnitude. Note that the aperture correction is negligible 
for our aperture sizep] 

We also obtained some exploratory exposures during an ob- 
servi ng run at the NTT using the tri-band ULTRACAM im- 
ager ( Dhillon et al. 2007 1 and managed to obtain an additional 
z-band and g-band image for PSR J2256-1024. The data pro- 
cessing was performed using the ULTRACAM data reduction 
pipelin^*^ 

3. RESULTS 

We found optical counterparts to all four pulsar binaries 
(see Figure [T]). The association of the optical counterparts 
with the irradiated pulsar companions was confirmed in all 
cases by variability at the known orbital periods, as can be 
seen in the light curves shown in Figure |2] (see also online 
Tables |2]|3]|4]and|5]l. 

As discussed in Section [T| the light curves can be used to 
constrain several parameters of the system. The orbital mod- 
ulation is caused by a combination of irradiation, which pro- 
duces a maximum at the superior conjunction of the com- 
panion, and ellipsoidal variations, which yield maxima at 
quadrature. From the colors at superior and inferior con- 
junction of the companion, one can constrain the dayside and 
nightside temperatures and thus irradiation. For poorly sam- 
pled data like ours, the nightside temperature cannot be con- 
strained very well, since it is degenerate with the orbital in- 
clination. Good coverage is required to distinguish between 
the high-inclination case of seeing only the brighter, unirradi- 
ated nightside at inferior conjunction of the companion, or the 
low-inclination case where one observes a dimmer nightside 
combined with a sliver of the irradiated side. However, for 
strongly irradiated systems, where the dayside is much hotter 
than the nightside, the irradiated temperature can be deter- 
mined relatively securely from just the colors near superior 
conjunction of the companion. The amplitude of ellipsoidal 
variations directly depends on the filling factor of the com- 
panion. Since the amplitude is small relative to irradiation, 
one must possess good quality sampling near quadrature in 
order to constrain the filling factor accurately. However, the 
physical size of the star is also directly related to its apparent 
luminosity given its distance, and thus any independent in- 
formation on the latter, such as from the dispersion measure, 
can be used to constrain the filling factor For small filling 
factorsp^below ^0.55, the scaling between filling factor and 
distance is linear (as expected), but for larger filling factors 
its dependence becomes sub-linear, as the ellipsoidal varia- 
tions and gravity darkening start to counteract the irradiation 
at maximum light (superior conjunction), leading to a ^ 20% 
reduction relative to the linear extrapolation when the com- 
panion fills its Roche lobe. Finally, we note that constraining 
parameters such as inclination and filling factor is generally 
limited by the quality of the sampling and the uncertainty on 
individual data points, whereas the determination of the tem- 
peratures are mostly affected by the absolute calibration of 
each band. 

*^ |ht tp ://heasarc.gsfc.nas a. gov/docs/swif t /| 

|analysis/uvo"t_ diqest / apercor . html 

^" |http : //deneb .astro . Warwick ■ ac .uk/phsaap/l 

|software/| 

' ' We define the filHng factor as the ratio of the volume-averaged stellar 
radius and the volume-averaged Roche-lobe radius. 



To get first-order estimates of the system parameters, we 
modelled the li ght curves of these irradiated binaries us- 
ing I car u;^ ( Breton et al."201l') along with BTSettl at- 
mosphere models (Allard et al. 2003 2007 2010) available 
through the Phoenix web simulato^^ The free parameters in 
the model are the mass ratio q, the orbital inclination /, the 
distance d and reddening Ay to the system, the companion's 
filling factor f = R/R]^ (where R is the volume-averaged ra- 
dius and /?L the volume-averaged radius of the Roche lobe), 
and its dayside and nightside temperatures Tday and T^^ight. We 
assumed that the companion is tidally locked to the pulsar 
and added Gaussian priors on the distance and reddening. We 
set the mean of the distance priors to the valu es inferred from 
dispersion measure using the NE2001 model (Cordes & Lazio 
2002| l, and, since our sources are well outside the dust layer. 



we estimated the reddening as be ing equal to the total along 
each line of sight as inferred by [Schlegel et al.| (11998) (see 
values in Table [TJ. In both cases, we took standard devia- 
tions corresponding to 25% relative uncertainties. Also, since 
the mass ratio is unconstrained without spectroscopic infor- 
mation, we imposed that the mass of the pulsar lies in the 
range 1-3 Mq . For the model fitting, we used the uncertain- 
ties from Tables [2] [3] |4] and [5] on individual data point, and 
a band-to-band uncertamty corresponding to our band cali- 
bration errors to which we added in quadrature an additional 
0. 1 mag to account for systematics in the zero-point calcula- 
tion inherent to the atmosphere models (this is important only 
for PSR J2256-1024 for which the photometry could be accu- 
rately tied to SDSS). In the following subsections, we discuss 
the main results of the light curve modelhng (see Figure[2]for 
best-fit light curves). 

3.1. PSR J0023+0923 

PSR J0023H-0923 stands out from other black widows in its 
relative long spin period, and thus relatively low spin-down 
luminosity, which makes the expected irradiation effects on 
the companion relatively mild. The data constrain the col- 
ors at maximum and minimum light well, with g-i^ 1 .5 and 
3.0, respectively. These correspond to flux-averaged tempera- 
tures of ^4500 and 3200 K, respectively, close to what we in- 
fer from our models: 4800 ± 2000 and 2900 ± 700 K, respec- 
tively. Note that the latter dayside temperature is strictly for 
the part directly facing the pulsar, hence quite a bit higher than 
the flux-averaged temperature inferred from the color. The 
substantial change in color shows that much of the irradiated 
side must be hidden at inferior conjunction of the companion, 
and our models demonstrate that this excludes orbital inclina- 
tions lower than / < 40°; the preferred value is / = 58 ± 14°. 

The filling factor appears relatively small at 0.30 ± 0.30, 
though one must be cautious since the probability distribution 
is highly non-Gaussian with a mode at ^ 0. 15 and an extended 
heavy tail all the way to a filling factor of unity (median value 
at 0.40). This could imply a much smaller companion than in- 
ferred for the prototype black widow, PSR B1957-H20, which 
has a filling factor fairly close to unity ( [Reynolds et al.|2007] ). 
Indeed, the nominal implied size, of ^O.O57?0 is smaller 
than is possible for a solar-composition object, suggesting 
that the filling factor is on the upper end of our inferred range 
(and thus the distance somewhat larger than inferred from the 
dispersion measure; see Section |4]i. However, indirect evi- 
dence that the companion is not close to filUng its Roche lobe 
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Figure 1. Finding cliart of the four irradiated pulsar companions presented in tliis paper, observed in i'-band with GMOS-North. In each panel, the location of 
the counterpart is indicated with crosshairs. 

comes from the fact that this source shows no radio eclipses 
at 350 MH2p"] unlike most black widows (Hessels et al., in 
prep.; though we could also be seeing the system relatively 
face -on). 



likely though given that our lightcurve fitting yield similar or- 
bital incUnation ranges (see Table[T]i. 



3.2. PSR J2256-1024 

This system appears to be quite similar to PSR J0023-I-0923, 
showing slightly redder colors. Near maximum, the fortu- 
itously phased ULTRACAM points yield g-z — 2.5, while 
near minimum, the Gemini data yield g-i ~ 3.5, correspond- 
ing to flux-averaged temperatures of ^3700 and 3100K, re- 
spectively. Since for PSR J2256-1024 our photometric cal- 
ibration is more secure, our fits yield well-constrained day 
and nightside temperatures, of 4200 ± 700 and 2450 ± 350 K, 
respectively. The modelling favors an intermediate orbital in- 
clination, / = 68± 11°. 

Like for PSR J0023H-0923, the best-fit filUng factor is 
small, 0.40 ± 0.20, although since the orbit is somewhat 
wider, the inferred radius of 0.09 is not inconsistent 
with a solar-composition, degenerate object. Contrary to 
PSR J0023-I-0923, however, this system shows radio eclipses, 
suggesting that, perhaps, the size is underestimated because 
the distance derived from the dispersion measure is too small 
(see also Section |4|, and/or that it is simply closer to being 
observed edge-on. The latter explanation appears slightly un- 

Radio eclipses are strongly radio -frequency dependent, being more dra- 
matic at lower frequencies (see, e.g. [Archibald et al.|2009| [Stappers et al.| 
[T996) . 



3.3. PSRJ1810+1744 

Like PSRs J0023H-0923 and J2256-1024, this system is a 
canonical black widow. One should expect a larger irradia- 
tion than in the case of PSRs J0023H-0923 and J2256-1024, 
given the combination of a compact 3.6-hr orbit and a more 
energetic pulsar (in fact, this is the fourth fastest-spinning pul- 
sar known in the Galactic field). Its light curves are quite 
puzzling since we did not find a combination of parameters 
that successfully account for all the data points. Indeed, the 
/-band light curve in particular seems inconsistent with be- 
ing symmetric around minimum light, being dimmer before 
minimum than after Asymmet ric hght curves hav e also been 
seen for PSR J205 1-0827 (Stappe rs et "dl]|200l] l, though in 
that case the source was brighter before minimum than it was 
afterwards, and the light curve showed variations between dif- 
ferent sets of observa tions (in contras t to what is seen for 

[2007) l. In the g-band 



PSR B1957H-20; e.g., Reynolds et al 



light curve, the data point near minimum is significantly off- 
set from the predicted value. As we mentioned in Section [2j 
particular care was taken for the flux extraction of this source. 
Hence we believe that the pulsar companion was intrinsically 
faint at this orbital phase. We note that some irradiated pulsar 



companions, such as PSR J13 11-3430 (|Ro mani et al.l2012[ ), 
display flares so it is not impossible that PSR J 1 8 1 Oh- 1 744 
also shows variability. Further observations with better or- 
bital coverage should help understand the odd behavior of this 
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Figure 2. Light curves of tlie four irradiated pulsar companions. GMOS 
GMOS g, ULTRACAM z, ULTRACAM ^> and UVOT uvwl bands are marlced 
by red circles, green triangles, magenta plus-signs, green cross-signs, and 
blue diamonds, respectively. The black lines display the best-fit light curve 
calculated using Icarus. 

source. 

As a result of the poor fit, we can only give qualitative con- 
straints on the parameters. Given the large modulation, the 
inclination has to be relatively large, though the lack of flat- 
tening near minimum implies that the dayside has to remain 
partly visible; we find / ~ 48 ± 7°. Our best-fit model also 
yields a very hot dayside, with ^ 10000 K, and a much colder 
nightside, with ^3100±850K. Given the larger distance to 
the system, ^ 2 kpc, a larger filling factor can be accommo- 
dated. Our data are consistent with a Roche-lobe filling star, 
though the poor fit and the lack of photometric coverage near 
maximum yields large uncertainties, hence explaining a heavy 
tail in the distribution of possible filling factor at lower val- 
ues. The inferred filling factor of 0.80 ±0.30, which implies 
a companion radius of O.15/?0. 

Given the poor fit, however, we caution that the above re- 
sults cannot be taken at face value. In particular, while our 
data cover minimum and thus directly constrain the nightside 
temperature (for g-/ ~ 1.3, one infers a flux-averaged tem- 
perature of ^4600K), maximum is not covered and hence 
we have no direct constraints on the dayside temperature. 
Nevertheless, the dayside must be relatively hot, since the 
source has / ~ near quadrature, implying a flux-averaged 
temperature of ^7500K. For many black widows, the color 



does not vary strongly between quadrature and maximum - 
see PSRs J0023 +0923 and J2256-H1 024 above, as well as 
PSR B1957H-20 ( jCaUanan et al.|1995| l. Given a dayside color 
of ^ - / < 0, and taking into account that the dayside temper- 
ature is somewhat hotter than the flux-averaged one, we infer 
a lower limit to the dayside temperature of > 8000 K. We will 
use this limit below. Using these values also reduces the irra- 
diation efficiency from the formal best-fit value of 1 .5 to 0.60, 
which make it more plausible energetically. 

3.4. PSR J2215+5135 

Unlike the systems we have already discussed, 
PSR J22 15-1-5 135 has a significantly more massive companion 
and is best characterized as a redback, with orbital properties 
similar to the prototype, PSR J1023H-0038 ( [Archibald et al. 
2009] l. Like PSR J1023H-0038, it shows only modest bright- 
ness variability. The Gemini data give only a single color, 
g-i ~ 0.3, implying a flux-averaged temperature of ^6600K 
near quadrature. The presence of UVOT data proves to be 
valuable at constraining parameters of this system, with the 
caveat that ultraviolet data are sometimes contaminated by 
chromospheric emission if the companion is magnetically 
active. The system appears to be viewed at an intermediate 
orbital inclination 70°) and a dayside temperature of 
6200 ± 500 K. Like other redback companions, the nightside 
is much hotter (^4800K) than that of black widow compan- 
ions. The filling factor is very tightly constrained to being 
Roche-lobe filling, which is consistent with the occurrence of 
radio eclipses for ^ 50% of the orbit (Hessels et al., in prep). 

4. DISCUSSION 

As can be seen in Figure [2] the model light curves gen- 
erally agree well with our data. While our light curves are 
not sufficiently well-sampled to produce strong constraints on 
the system parameters, they allow us to address two interest- 
ing aspects: the efficiency of the irradiation and the extent to 
which the companion fills its Roche lobe. 

To estimate the efficiency of the irradiation, i.e., the ef- 
fective fraction e\„ of the spin-down luminosity incident on 
the companion that is absorbed and re-radiated, we assume 
that the irradiating flux is thermalized and reradiated lo- 
cally, i.e., that it simply adds to the intrinsic flux wherever 
it impinges. Then, from the hottest point, one can derive 
a characteristic "irradiation temperature" = — r^^jg^t, 
which is related to the pulsar's spin-down luminosity, L^^, as 
ein-i-sd = 47rfl^(T7jJ*j., where a is the orbital separation and a is 
the Stefan-Boltzmann constant. Note that we implicitly as- 
sume that the pulsar spin-down energy is isotropically radi- 
ated. 

In Table [T] we list our inferred irradiation temperatures as 
well as the implied irradiation efficiencies. With the excep- 
tion of PSR J1810-I-1744, which we will discuss below, we 
find that the typical values of the irradiation efficiency of 
the new systems presented in this paper are consistent with 
those of PSRs B1957H-20 and J1023H-0038. Moreover, the 
other known irradiated p ulsar systems (PSRs J205 1-0827, 
J131 1-3430, J2339-0533[Stappers et al. '2001 ; Roman i et al.| 



|2012[ Romani & Shaw|201 \ \ also display similar efficiencies 
(^^%, ~30% and ~ 15%, respectively^] We conclude that 
the typical irradiation efficiency factor in these systems lies in 

^' PSR J205 1—0827 has a somewhat large value efficiency. The efficiency 
of PSR J2339— 0533 is based on the radio timing from Ray et al. (in prep.). 
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the range 10-30%, with 15% being a representative figure. 
Note that the intrinsic spread in values might be smaller if one 
compensates for the fact that the orbital separation is poorly 
known given the uncertainties in the inclination, the pulsar 
mass and its moment inertia. 

The energetics derived from the irradiation of the compan- 
ions are consistent with the idea that the relativistic wind, 
which is powered by the rotational spin down of the neutron 
star, is the major driver of the heating mechanism. The case 
of PSR J18 10+1744 is, however, puzzling. Its rather large in- 
ferred irradiation efficiency not only departs from the other 
known irradiated pulsar systems but it also implies an input 
energy larger than the nominal spin-down luminosity (though 
our analytic estimate of the lower limit of ~ 0.60 is below 
unity). One cause may be that our assumption of an isotropic 
wind is not justified. For instance, if the pulsar is aligned with 
the orbit and emits its wind preferentially in the equatorial 
plane (as is the case for, e.g., the Crab and Vela pulsars), the ir- 
radiation efficiency would be reduced. From our Gemini data, 
there is a clear indication that both g and / band light curves 
are not symmetric, with the companion being brighter after its 
inferior conjunction (phases ^ 0.0-0.15) than before (phases 
^ 0.85- 1.0). This could be an indication of non-isotropic 
heating or heat redistribution at the surface of the star. More 
detailed light curves of this system would help resolve these 
issues. 

Radio eclipses are observed in three out of four of these 
systems (Hessels et al. in prep.; Stairs et al. in prep.) and 
large increases in the dispersion measure at the ingress and 
egress in similar systems indicates total intra- binary electron 



column densities of A^^, ^ lO'^cm ^ (see, e.g.. 



Fruchter et al. 



1988b|). Given that eclipses are coincident with the inferior 



conjunction of the companions, plasma must certainly be sur- 
rounding them and hence some form of mass loss from the 
companion is required. In the case of a nearly Roche-lobe 
filling star, material is loosely bound to the surface and can be 
peeled off easily when exposed to a relativistic pulsar wind. 
Our work suggests that, however, some of the irradiated pul- 
sar companions that we have studied (PSRs J0023+0923 and 
J2256-1024) are not close to filling their Roche lobe. 

As mentioned in the previous section, a first possibility is 
that the distances inferred from the dispersion measures of 
these pulsars are underestimated. It is now well established 
that for pulsars located far off the Galactic plane, the mea- 
sured parallactic distances tend to be larger than the DM dis- 
tances (see, e.g., Gaensler et al. 2008 ; Roberts 2011). Roberts] 
( [201 1) shows that based on 13 sources with measured parallax 
and Galactic latitude larger than 10°, t/oivi/'^paraiiax = 0.66 ± 
0.26. Since the filling factor is correlated with the distance, 
we ran another set of fits using priors on the distance rescaled 
using the above conversion factor and error As a result, we 
found that the typical filling factor for PSRs 10023+0923^ 
J1810+1744 and J2256-1024 did not change significantl>0 
This comes from the fact that the dayside temperature are not 
very precisely constrained due to the large systematic uncer- 
tainties in the absolute calibration of the bands. Consequently, 
the larger distance priors tends to increase the dayside tem- 
peratures rather than changing the filling factors as one would 
expect. We also ran another set of fits, this time by holding 
the filling factor of the companions to unity and removing the 
distance priors in order to see how much further these sys- 
tems would need to be located in order to match the observed 

PSR J2215+5135 is already constrained to be Roche-lobe filling. 



fluxes. We found that the DM distances would need to be 
off by a factor 8.5, 1.2 and 2.3 for the three above sources, 
respectively. 

While it is not excluded that our distance estimates 
are wrong, it appears unlikely that PSRs J0023+0923 and 
J2256-1024's companions are Roche-lobe filling. Whether 
nearly Roche-lobe filling stars are required in order to explain 
the radio eclipses of the pulsars is uncertain since neither the 
mechanism supplying particles to the plasma nor the role of 
the pulsar at triggering it are understood. Caution should be 
taken before drawing definitive conclusions and precise dis- 
tance constraints from parallactic measurements would help 
shed light on this. Better-quality multi-color light curves will 
improve the measurement of the orbital inclination and ad- 
dress the contribution of ellipsoidal variations, hence also help 
to constrain the filling factor. It is worth mentioning that the 
second black widow system to be found, PSR J205 1-0827, 
also displays puzzling light curves. Previous work highlighted 
an ambiguous behavio r either indicating a fi lUng factor near 
unity or closer to 50% ("St appers et al.|200T) . 

A common feature of black widow and redback systems is 
the presence of non-secul ar orbital period derivatives in the ra- 
dio timing (see, e.g. Arzo umanian et al.|1994[|Lazaridis et al.| 
201 1 , Archibald et al. 200^. It has been suggested that grav- 
itational quadrupole couplmg of the compa nion with the or- 
bit might explain the orbital variability ( Applegate & Shaham 
il994j). The dissipated tidal energy would drive convection, 
which would power a dynamo-induced magnetic field and 
provide a significant source luminosity. While such a mech- 
anism would be compatible with the amplitude of the orbital 
variatio ns of PSR B 1957+20 and the luminosity of its com- 
panion ( Applegate & Shaham||1994] l, it would require some 
fine tuning - namely a ^50% filling factor - i n order to work 
for PSR J205 1-0827 (Lazari dis et aL][20n] l. These newly 
discovered irradiated pulsar systems could therefore provide 
extended leverage to test the gravitational quadrupole theory, 
since the Roche-lobe under-filling systems are predicted to 
display smaller orbital variability. 

5. CONCLUSION 

Our view of the binary pulsar population is currently shift- 
ing toward a new paradigm. Until the launch of Fermi, the 
bulk (^ 90%) of the known population i n the Galactic field 
consisted of pulsar-white dwarf systems (Lorimer & Kramer 
2004| l, while the remaining pulsar binaries had neutron star, 
mam sequence, very low-mass star or planet companions. 
Only about four of the known binary pulsars in the Galactic 
field were irradiated systems like those presented here. As of 
today, the number of irradiated systems has increased to over 
twenty members and candidates, which implies that they now 
account for about 10% of the binary pulsar populat ion out- 
side of g lobular clusters (based on the ATNF catalog, |Manch-j 
|ester et al.|20 05 ). It is clear that a large selection bias against 
finding these binaries in classical radio surveys existed until 
high-energy missions were added in the picture — and yet 
they still remain challenging to find. New radio pulsar sur- 
veys, benefiting from multibeam receivers and larger band- 
width, are also contributing to finding irradiated systems in 
blind searches, as was the case for PSR J2256-1024. Black 
widows and redbacks therefore constitute a fundamental com- 
ponent of the pulsar ecosystem, a component th at dominates 
among the fastest-spinning MSPs ( Hessels|2008 1 



The work presented here shows that the spin-down lumi- 
nosity of pulsars is a good indicator of the level of irradiation 
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sustained by their companions. We found that these systems 
display a rather universal irradiation efficiency ~ 10-30% 
for reprocessing the incoming energy flux from the pulsar's 
spin-down into heat on the companion's surface. As a result, 
one may easily estimate the brightness and amplitude of the 
optical light curves due to irradiation in these pulsar binaries, 
provided an orbital separation, if the pulsar spin-down lumi- 
nosity is known from timing or, alternatively, from the 7-ray 
luminosity. 

The typical reflection albedo of stars with temperatures in 
the range 2000- 10000 K is between 0.5 and 1.0 (Claret 2001) 
for atmospheres that are convective and in radiative equi- 
librium, respectively. Given the above 15% irradiation ef- 
ficiency, the above albedos imply that 10-30% of the en- 
ergy from the spin-down luminosity would actually reach 
the companion. It is worth noting that the temperature of 
PSR J1810H-1744 suggests that its outer envelope might be 
radiative, as opposed to convective in the other systems pre- 
sented here. If so, the observed irradiation efficiency would be 
a factor 2 larger because of the difference in the stellar albedo 
and this could partly explain why it appears unusual. If fur- 
ther studies of these systems find bow shock nebulae (like for 
PSR B1957H-20), it would allow an independent measurement 
of the energy loss by the pulsar, which would make for an 
interesting comparis on with that inferre d from the irradiated 
companion (see, e.g., van Kerkwijk & Ingle 2008 ). 

The possibility that some of these pulsar companions do 
not fill their Roche lobe leads one to ponder the underlying 
cause for radio eclipses in these systems. What is the mech- 
anism responsible for replenishing the plasma responsible for 
the eclipses? Other missing pieces of the evolutionary puzzle 
are: Were these companion stars closer to filling their Roche 
lobe in the past or have they contracted thermally since the 
mass transfer episode has terminated? Can black widows and 
redbacks completely destroy their companion and eventually 
become the isolated millisecond pulsars that we observe? Or, 
are we simply observing the systems that are incapable of de- 
stroying their companions on a relatively short timescale? 

The fact that some of these irradiated systems are relatively 
bright also offers an interesting opportunity for spectroscopic 
follow-up. Radial velocity curves should allow for the mea- 
surement of the component masses, and test whether these 
neutron stars typically are more massive, as found for the sin- 
gle system studied so far Further, detailed photometry work 
will also allow one to investigate the possible asymmetry in 
the Hght cui-ves of PSR J1810H-1744 and the reason for its 
anomalously large dayside temperature. 
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Table 2 

Photometry of PSR J0023+09 



Time 


Orbital Phase" 


Hux 


Hux Error'^ 


Magnitude'' 


Magnitude Error^ 


(MJD) 




(hJy) 


(hJy) 






i'-banU, CMOS 


55445.513890 


0.6410 


0.989 


0.065 


23.908 


0.070 


55449.575790 


0.9056 


0.024 


0.057 


26.756 


0.511 


55449.587001 


0.9863 


0.001 


0.065 


26.964 


0.599 


55449.612403 


0.1693 


0.106 


0.062 


26.099 


0.428 


i-band, GMOS 


55445.507793 


0.5970 


6.397 


0.252 


21.883 


0.043 


55445.519981 


0.6848 


4.368 


0.224 


22.296 


0.055 


55449.569694 


0.8616 


0.956 


0.239 


23.889 


0.244 


55449.581881 


0.9494 


0.637 


0.240 


24.266 


0.330 


55449.606307 


0.1254 


1.715 


0.238 


23.295 


0.146 


55449.618493 


0.2132 


3.027 


0.242 


22.691 


0.086 



"Orbital phases are measured from the companion's inferior conjunction. 

AB magnitudes in the Lupton system, m = njQ — 2.5 log fo' - (2.5 log e) sinh"! (f/2b'), using softening parameters b' = 0.059 and 0.233 /i7y for the g and band, 
respectively, and zero-point = —48.6. 

'^The flux and magnitude errors represent the formal uncertainties. One should also add in quadrature the zero-point caUbration errors, which are 0.319 and 
0.757 mag in g and i band, respectively. 



Table 3 

Photometry of PSR J2256-10 



Time 


Orbital Phase" 


Hux 


Flux Error^ 


Magnitude'' 


Magnitude Error^ 


(MJD) 




(iJjy) 


(iJjy) 






g-band, GMOS 


55428.546679 


0.2850 


1.871 


0.060 


23.219 


0.035 


55449.283273 


0.6936 


1.835 


0.029 


23.241 


0.017 


55454.418742 


0.8170 


0.284 


0.026 


25.256 


0.098 


55469.350249 


0.9567 


0.020 


0.024 


27.478 


0.483 


i-band, GMOS 


55428.540594 


0.2564 


6.770 


0.130 


21.823 


0.021 


55428.552779 


0.3136 


10.368 


0.125 


21.360 


0.013 


55449.277175 


0.6649 


9.902 


0.088 


21.410 


0.010 


55449.291633 


0.7328 


5.578 


0.069 


22.033 


0.013 


55454.412643 


0.7884 


3.321 


0.078 


22.596 


0.026 


55454.424831 


0.8456 


1.559 


0.074 


23.416 


0.051 


55469.344150 


0.9280 


0.867 


0.067 


24.049 


0.083 


55469.356340 


0.9853 


0.816 


0.062 


24.114 


0.081 


g-band, ULTRACAM 


55323.421265 


0.4664 


3.979 


1.606 


22.245 


0.335 


z-band, ULTRACAM 


55323.421265 


0.4664 


45.700 


5.609 


19.733 


0.129 



"Orbital phases are measured from the companion's inferior conjunction. 

''AB magnitudes in the Lupton system, m = — 2.51og^)' — (2.51oge)sinh"'(//26')j using softening parameters b' = 0.025, 0.064, 1.673 and 5.845 ^iJy for the 
g (GMOS), i, g (ULTRACAM) and z band, respectively, and zero-point m„ = --48.6. 

'^The flux and magnitude errors represent the formal uncertainties. One should also add in quadrature the zero-point caUbration errors, which are 0.104, 0.119, 
0.010 and 0.010 mag in g (GMOS), i, g (ULTRACAM) and z band, respectively. 
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Table 4 

Photometry of PSR J1810+17 



Time 


Orbital Phase^ 




Hux Error^ 


Magnitude'' 


Magnitude Error^ 


(MJD) 




(/<./v) 


(/i,/v) 






g-band, CMOS 


55442.236873 


0.7092 


40.430 


0.119 


19.883 


0.003 


55443.292770 


0.8355 


11.639 


0.053 


21.235 


0.005 


55446.293441 


0.0870 


6.537 


0.061 


21.861 


0.010 


55447.310477 


0.9509 


1.066 


0.049 


23.828 


0.050 


i-band, GMOS 


55442.230765 


0.6680 


42.101 


0.141 


19.839 


0.004 


55442.242953 


0.7502 


29.091 


0.124 


20.240 


0.005 


55443.286671 


0.7943 


20.585 


0.132 


20.616 


0.007 


55443.298861 


0.8766 


8.369 


0.107 


21.593 


0.014 


55446.287344 


0.0458 


5.916 


0.120 


21.969 


0.022 


55446.299530 


0.1281 


14.863 


0.126 


20.969 


0.009 


55447.304378 


0.9098 


5.780 


0.095 


21.995 


0.018 


55447.316567 


0.9920 


3.715 


0.107 


22.474 


0.031 



"Orbital phases are measured from the companion's inferior conjunction. 

''AB magnitudes in the Lupton system, m = -2.51ogfo'-(2.51oge)sinh~'(//2fe'), using softening parameters b' = 0.051 and 0.099 ijjy for the g and band, 
respectively, and zero-point nij, = —48.6. 

'^The flux and magnitude errors represent the formal uncertainties. One should also add in quadrature the zero-point caUbration errors, which are 0.301 and 
0.308 mag in g and i band, respectively. 



Table 5 

Photometry of PSR J2215-H51 



Time 


Orbital Phase'' 


Hux 


Flux Error' 


Magnitude* 


Magnitude Error' 


(MJD) 












g-band, GMOS 


55531.246057 


0.1970 


82.000 


0.232 


19.115 


0.003 


i-band, GMOS 


55531.192316 


0.8855 


59.559 


1.180 


19.462 


0.022 


55531.238446 


0.1529 


71.956 


0.607 


19.257 


0.009 


55531.252146 


0.2323 


95.277 


0.448 


18.952 


0.005 


55531.256359 


0.2567 


102.023 


0.443 


18.878 


0.005 


uvwl-band, UVOT 


55399.269193 


0.1227 


0.670 


2.759 


22.979 


0.749 


55399.335860 


0.5092 


4.868 


2.407 


22.036 


0.417 


55399.402162 


0.8936 


-1.056 


2.269 


23.450 


0.603 


55399.536221 


0.6707 


2.450 


2.136 


22.525 


0.499 


55399.602885 


0.0572 


-1.364 


2.157 


23.531 


0.561 


55399.665753 


0.4216 


5.503 


2.026 


21.930 


0.325 


55399.730436 


0.7966 


2.269 


2.144 


22.568 


0.511 


55399.799357 


0.1961 


-0.123 


1.893 


23.196 


0.521 


55399.865723 


0.5808 


5.094 


2.032 


21.997 


0.342 


55399.933775 


0.9753 


-0.077 


1.976 


23.183 


0.544 



^Orbital phases are measured from the companion's inferior conjunction. 

''AB magnitudes in the Lupton system, m = - 2.5 log b' - (2.5 log e) sinh"' {f/2b'), using softening parameters b' = 0.242, 0.462 and 1.972 fUy for the g, i 
and uvwl band, respectively, and zero-point = -48.6. 

'The flux and magnitude errors represent the formal uncertainties. One should also add in quadrature the zero-point caUbration errors, which are 0.302 and 
0.302 mag in g and band, respectively. 



